hep-th/0304079 
PUPT-2080 



Proton Decay in Intersecting D-brane 

Models 

Igor R. Klebanov 1 and Edward Witten 

Institute for Advanced Study, Princeton, NJ 08540, USA 



Abstract 

We analyze proton decay via dimension six operators in certain GUT-like models 
derived from Type IIA orientifolds with D6-branes. The amplitude is parametrically 

— 1/3 

enhanced by a factor of ot GJ j T relative to the coresponding result in four- dimensional 
GUT's. Nonetheless, even assuming a plausible enhancement from the threshold cor- 
rections, we find little overall enhancement of the proton decay rate from dimension 
six operators, so that the predicted lifetime from this mechanism remains close to 10 36 
years. 
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1 Introduction 

Grand Unified Theories in four dimensions have had impressive successes, in accounting 
for the quantum numbers of fermions [1], in predicting - with the aid of supersymmetry 
- a value of sin 2 6 W that is in excellent agreement with experiment [2] , and in pointing 
twenty years in advance of decisive measurements to the right order of magnitude of 
neutrino masses [3]. GUT theories also make the exciting prediction that the proton 
may decay with a lifetime close to present experimental bounds. 

With supersymmetry, even if one assumes an i?-parity symmetry to avoid catas- 
trophic proton decay at low energies, there actually are two different GUT-based mech- 
anisms for proton decay. The proton may decay due to dimension five operators of the 
form / d 2 9Q 3 L, where Q and L are quark and lepton superfields. Or it may decay by 
dimension six operators of the form / d 4 9Q 2 Q*L* coming from gauge boson exchange; 
this effect is similar to proton decay in GUT's without supersymmetry, but is signif- 
icantly slower because supersymmetry raises the GUT scale. In the simplest models, 
proton decay by dimension five operators dominates, and in fact present experimental 
bounds make life difficult for these models [4, 5, 6] (see, however, [7]). This suggests 
the possibility that some mechanism suppresses the dimension five operators and that 
gauge boson exchange may after all be the dominant mechanism. Various methods 
to suppress the dimension five operators while preserving the successes of GUT's are 
known, though they tend to be somewhat technical. For example, one construction 
based on discrete symmetries [8] is fairly natural in the class of models that we will 
consider in this paper. 

Most GUT-like string-based models of particle physics do not precisely lead to four- 
dimensional GUT's, since unification takes place only in higher dimensions (for an early 
review, see [9], chapters 15 and 16), leading among other things to possibilities for GUT 
symmetry breaking by discrete Wilson lines [10] and to a higher-dimensional mechanism 
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for doublet-triplet splitting [11]. Generally, models such as the heterotic string on a 
Calabi-Yau manifold lead to qualitatively similar issues concerning proton decay to 
those in GUT's, though the details are somewhat different. One important reason that 
the details are different is that because of the higher-dimensional structure, there is 
generally in these models no precise answer to the question, "Which color triplet gauge 
boson is the SU (5) partner of the standard model gauge bosons?" There is a lightest 
color triplet gauge boson, whose wave-function in the compact dimensions depends on 
the details of the model, and its exchange may, depending on the model, give the right 
order of magnitude though not the correct numerical value of the dimension six part 
of the proton decay amplitude. 

Recently [12], these issues were reconsidered in another class of models - M-theory 
on a G 2 manifold. (Actually, in many cases these models give dual descriptions, useful 
in a different region of the parameter space, of the same models that can be studied 
via the heterotic string on a Calabi-Yau manifold.) The GUT threshold correction to 
proton decay was computed in this class of models and was seen to give a potential 
enhancement of the proton decay rate. It was also shown that in this type of model, be- 
cause of the way quarks and leptons are localized, exchange of the lightest color triplet 
gauge boson does not dominate the proton decay amplitude. On the contrary, a field 
theoretic attempt to compute a proton decay amplitude by summing over Kaluza-Klein 
harmonics runs into an ultraviolet divergence, and consequently the correct answer de- 
pends on the cutoff provided by M-theory. Formally, this UV divergence enhances 

— 1/3 

the proton decay amplitude by a factor of a GL } T compared to what it would be in a 
four-dimensional model. (This factor must be combined with the threshold factor, of 
course.) The coefficient of ct^j^ is universal and independent of the details of the 
model (such as how the GUT symmetry is broken) - in fact, the dominant four-fermi 
operator is local and invariant under the full GUT symmetry, in contrast to the usual 
situation in four-dimensional GUT's. However, present knowledge of M-theory does 
not make it possible to compute the numerical coefficient of this operator. 

As an alternative, we will in the present paper consider another dual class of mod- 
els in which the calculation can be done. These models fall in the general class of 
intersecting brane worlds [13, 14, 15, 16, 17, 18, 19] - Type IIA superstrings with 
gauge bosons supported on D6-branes and chiral matter multiplets at intersections of 
the D6-branes. Many supersymmetric GUT-like models have been constructed along 
these lines [20, 21, 22], though so far none with the mechanism of symmetry breaking 
by Wilson lines on the _D6-brane that we will assume in the present paper to justify 
borrowing various above-cited results. 1 (In actual model-building, one typically con- 

1 In the dual context of M-theory on G 2 manifolds, some examples of symmetry breaking by Wilson 
lines were considered in [23]. 
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siders branes in an orbifold of a six-torus, but more generally one may start with any 
Calabi-Yau manifold: see, for example, [24].) Threshold corrections in intersecting D6- 
brane GUT-like models were calculated in [25]. These models are dual to M-theory 
on a manifold of G2 holonomy as g s , the Type IIA string coupling constant, becomes 
large. The chiral matter fields are localized (at brane intersections) in a similar way to 
what happens in M-theory on a G2 manifold, and accordingly we will find the same 
anomalous factor of a Gl j T in the proton decay amplitude. The difference is that in 
perturbative Type IIA string theory, everything is explicitly calculable, and hence we 
will be able, at least for g s <C 1, to be precise about the numerical factors. Regrettably, 
the numerical factors accompanying a^j^ are such that even assuming a plausible en- 
hancement from the threshold factors, the proton decay rate is comparable to or only 
slightly greater than that in standard SUSY-GUT's. Of course, the precise factors that 
appear in the corresponding M-theory (or large g s ) limit are still unknown. 

In this paper we will not be concerned with any specific model, but will rather try 
to incorporate universal features of the GUT-like intersecting D-brane models relevant 
to calculation of the proton decay rate. Typically one needs a stack of D6-branes 
intersecting an orientifold fixed sixplane along 3 + 1 dimensions. On the covering 
space we then have a stack of D6-branes intersecting an image set of D6'-branes. If 
there are 5 D6-branes in the stack, then on the covering space we find gauge group 
SU(5) x £[/(5), with open strings localized at the intersection transforming in the 
bifundamental representations (5, 5) + (5, 5). After the orientifold projection we find 
SU(5) gauge theory with matter in 10 + 10. Our goal is to calculate the 4-fermion 
contact term whose £77(5) structure is 10 2 10 2 . Such an interaction mediates proton 
decay processes such as p — > 7r°e£ (and, depending on the assumed flavor structure, 
other modes with e\ replaced by \x\ and/or 7r° by K°). The calculation is conveniently 
carried out on the covering space where we need the 4-point function for two (5, 5) 
states and two (5, 5) states. The calculation is sensitive to the local structure of the 
intersection, and is insensitive to how the D6-branes are wrapped around the compact 
space, as long as its size is somewhat greater than the string scale so that worldsheet 
instantons are suppressed. 

Four-dimensional GUT's also have dimension six operators 10 10 5 5 which lead to 
proton decay modes such as p — > 7r°e^ and p — > tt + v; as observed in [12] in the 
analogous M-theory case, such interactions do not arise generically in a model of this 
type. To get such an interaction, we need two stacks of branes that are not mirror 
images to meet on an orientifold plane; in general, this would be a coincidence. Both 
four-dimensional GUT's and the brane worlds have additional 5 2 5 2 interactions (which 
in the brane worlds arise from brane intersections away from the orientifold planes), 
but these do not violate baryon number. 
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2 Vertex Operators 



Without a loss of generality we may assume that the D6-branes are oriented in the 
0123468 directions. The D6'-branes intersect them along the 0123 directions that will 
have the interpretation of a 3 + 1 dimensional 'intersecting brane world.' To specify its 
orientation in the six transverse directions, we define the complex coordinates 

z\ = x A + ix 5 , Z2 = x 6 + ix 7 , Z3 = x 8 + ix 9 . (1) 

In order to preserve Af = 1 supersymmetry in 3 + 1 dimensions, the rotation must act 
as an SU(3) matrix on the three complex coordinates [13]. Choosing the matrix to be 
diagonal, we see that the rotation that turns the D6-branes into the D6'-brane acts as 

z\ — > exp(7r^i)zi , z 2 — > exp(iri6 2 )z 2 , z 3 — > exp(7ri9 3 )z3 , (2) 

where 

0i + 9 2 + #3 = 2 mod 2Z . (3) 

We need to construct the vertex operators for the 6 — 6' and 6' — 6 open strings. 
Such operators, when inserted on the boundary of the disk, create discontinuities in the 
boundary conditions for the transverse coordinates and their fermionic partners [26]. A 
standard method for studying correlators on the upper half-plane is the doubling trick 
where the half-plane is replaced by the entire plane, but with only the holomorphic part 
of the field on it (for a review, see [27]). If we consider complex fields X 1 corresponding 
to the transverse coordinates z\ then their Laurent expansion around the insertion 
of the vertex operator has mode numbers shifted by Qi [35]. This is analogous to 
what happens near a twist field introduced in orbifold theories and we can use similar 
methods for calculating the correlation functions. 

Consider a bosonic twist field cr+ which creates a discontinuity in a complex coordi- 
nate X. Then we have the OPE [28, 29, 30] 

dX(z)a+(w) ~ (2 - wY" 1 ^ , dX(z)a+(w) ~ (z - w)- e r' + . (4) 

The dimension of a + is A CT = 6(1 — 9)/2, and we will assume that lies between and 
1. Now consider the complex fermion vp, which is a worldsheet superpartner of X. By 
worldsheet supersymmetry, its mode numbers are also shifted by 9. In the Ramond 
sector, cr + must be accompanied by the fermionic twist s + such that 

*/>(z)s+(w) ~(z- w) e - l/ H + , tf(z)s+(w) ~ (z - wf/ 2 - e t' + . (5) 

If we bosonize ip into a real scalar H, whose Green's function is 

{H(z)H(w)) = -log(z-w) , (6) 
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then 

i> = exp(iH) , i) = exp(-iH) , s + = exp[i(6 - 1/2)H] . (7) 
This means that the dimension of s + is 




Therefore, 

A s + A a = 1/8 , (9) 

independent of the rotation angle 9. 

The complete vertex operator that creates a massless 6 — 6' string in the R sector is 

3 

V+ = \ i I u a e-^ /2 S a e lk ^ II« S V) ( 10 ) 

where <fi is the bosonized ghost field. If we bosonize the fermions ip 11 , \i = 0,1,2,3, 
into scalars hi and I12, then the 3 + 1 dimensional spinor S a ~ exp[i(sihi + S2/12V2], 
where s±, s 2 = ±1. The GSO projection for 6 — 6' strings requires that sis 2 = 1, which 
restricts us to a spinor of definite chirality. A} is a Chan-Paton factor with an index 
i in the fundamental of the first SU(5), and an index / in the antifundamental of the 
second SU(5). Since /i(e~^ 2 ) = 3/8 and h(S a ) = 1/4, we see that the requirement 
that the dimension of V+ is 1 gives the 4-d massless dispersion relation k 2 = 0. 

The vertex operator V- for 6' — 6 strings is constructed similarly, except we need to 
replace the 0+ and s + twists by cr_ and s_. The latter are defined by sending 6 — > 1 — 6 
in the OPE. The GSO projection now requires that Sis 2 = — 1; this corresponds to 4-d 
spinor chirality opposite to that of the 6 — 6' string. Thus, we have the vertex operator 

V_ = \ I i u & e-*' 2 S d 'e ik » x ' 1 {[(atsi) . (11) 

3 The 4-fermion Amplitude in String Theory 

Our goal is to calculate the 4-fermion amplitude 

f 1 dx(V-(0)V + (x)V-(l)V + (oo)) . (12) 
Jo 

The crucial ingredient in this calculation is the knowledge of the correlator of the 
bosonic twist fields a. In the orbifold theories, these were calculated in [28, 29, 30]. 
We will take the square root of their result to account for the fact that we have the 
boundary twist fields rather than the bulk ones. (This intuitively plausible result has 
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been derived in [31] and, up to normalization, in [32]; see also [33, 34] for a treatment 
of cases where the intersecting D-branes have different dimensionalities.) This gives 



(<7_(0)<7 + (a;)<7_(l)<7 + (oo)) ~ 0rin(7r0)[a;(l - x)]~ 2A ° [F{x)F{\ - x)]^ 2 , (13) 
where 

F(x) = F(6,l-6;l;x) (14) 

is the hypergeometric function and ~ represents a numerical constant. 
On the other hand, the correlator of the fermionic twists is simply 

(s_(0)s + (x)s_(l)s+(oo)) ~ [x(l-x)]- 2A ° , (15) 

while the ghost factor contributes 

(e-^ 2 (0)e^ /2 (a;)e^ /2 (l)e"^ 2 (oo)) ~ [x(l - x)]- l ' A . (16) 

The correlator of the 4-d spin fields multiplied by spinor polarizations gives an x- 
independent factor Tii 7^2^37^ ^4- This x-independence is evident from the form of the 
bosonized correlator 

(exp[i(/ii - h 2 )/2] (0) exp[i(/n + h 2 )/2] (x) exp[-i(/n - h 2 )/2]{l) exp[-i(/n + h 2 )/2] (oo)) . 

(17) 

Putting together all the parts of the vertex operators, defining as usual 

s = -(h + k 2 ) 2 , t = -(k 2 + h) 2 , (18) 
and including both s- and t-channel diagrams, we find the amplitude 

4 

A A (ki, k 2 , k 3 , k 4 ) = iC(27r) 4 5 4 (^ h^^u^^^UA Tr (AiA 2 A 3 A 4 + AiA 4 A 3 A 2 ) 

rl 3 

• / dxx-^'^l - xy^'^JsmiTre^Fie.A - e f A-,x)F(6 l ,l - e t] i-,i - x)}- 1 / 2 , 

(19) 

where C is a normalization that we will fix later, while the A's and A's are Chan- 
Paton matrices of 6 — 6' and 6' — 6 strings. In four dimensions, with U\,u% having 
one chirality and u 2 , w 4 the other, we have ui^^u^^u^ = — mi7 m -u 4 M37 m -u 2 . This 
insures the antisymmetry of the amplitude under permutation of particles 2 and 4 (or 
of particles 1 and 3). 
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We note that the hypergeometric functions are crucial for the convergence of the 
integral at the endpoints of the integration region for vanishing momenta fcj. As x — > 0, 

F(x) -> 1 , F(l - x) -> ln(5/x) , (20) 

7T 

where 

In 5(0) = 2^(1) - ^(0) - ^(1 - 0) . (21) 

Therefore, even with s = t = 0, we find a convergent integral. Its behavior near 
x = is 

^3/2 /" ^[l n (i/ x )]-3/2 (22) 

JO X 

The physical reason for the absence of an IR divergence at vanishing momentum is 
the following. Even though there is a massless intermediate state from the untwisted 
sector, the special kinematics of this problem gives no IR divergence. For example, if 
the s-channel intermediate state is a massless 6' — 6' string, which is a gauge boson, then 
it can carry arbitrary momentum q along the directions of the _D6'-brane orthogonal 
to the intersection. These three components of momentum have to be integrated over: 

/ d 3 q [ dxx a 'f- a ' s - x = 7r 3 /Vr 3/2 / dxx- l - a ' s [\n{l/x)}-^ 2 , (23) 
J Jo Jo 

which converges near x = even for s = 0. Thus, going from the effective field theory 
to the form of the string integrand near x = 0, we find the replacement 



f vr 3 /V)- 1/2 / dxx' l ~ a ' s [\n{l/x)}^/ 2 . (24) 

J q 2 — s Jo 



This observation also allows us to normalize the 4-fermion amplitude. If we consider 
the s-channel factorization on the 6' — 6' gauge boson, which is the leading term as 
x — > 0, then we should find 

A 4 (k u k 2 ,k 3 ,k 4 )^-t —— £ — e , (25) 

J {Zir) ' k A — te 

where the 3-point function for emission of a gauge boson of Chan-Paton indices J, / is 
A T j{ki, k 2 , fc 3 ) = iK^- s {a'fl\2TxY5\Y. fc^i7^2 (AiA 2 ) J j . (26) 



i=i 



The index /i corresponds to the polarization of the intermediate vector state. Since 
the gauge coupling on D6-branes is given (according to [35], eqn. 13.3.25) by g 2 m = 
(27r) 4 ^(a') 3/2 , we have 

K 2 = (2tt) 4 . (27) 



7 



Performing the integral on the RHS of (25) with the replacement (24), we find 

2Tr(A 1 A 2 A 3 A 4 )^(a / )(27r)5 4 (E^)«i7^2M37^47r 3/2 I x' 1 ^' s [\n{l/ x)}^ 2 . (28) 

Equating this to the contribution to A 4 , given in eqn. (19), from the s-channel 6' — 6' 
massless intermediate state, 

i Tr(A 1 A 2 A 3 A 4 )C(27r) 4 5 4 (^ k^u^u^i^' 2 f ar^Ml/x)]" 3 / 2 > ( 29 ) 

i=i ^° 

we find that 

C = (2n)- 3 K 2 g s a' = 2ng s a' . (30) 
The low-energy limit of the 4-fermion amplitude is hence, 

4 

A 4 (ki, k 2 , fc 3 , fc 4 ) =i(27r# s c/)^(#i, 6> 2 , 6> 3 )(27r) 4 5 4 (^ h)ui^u 2 Uz r )^u 4 

i=i , . 

r _ (31) 
Tr (AiA 2 A 3 A 4 + AiA 4 A 3 A 2 J , 

where 

/•l //t' 3 

m,6 2 ,e 3 )= / IJ Vsin(7rgi)[F(^, 1 - g i; 1; x)F(^, 1 - g <; 1; 1 - x)}' 1 ' 2 . 

Jo x{i x) i=1 

(32) 

This has been derived for intersections of two stacks of branes without any orientifold- 
ing. Thus, in an intersection of a stack of N D6-branes with a stack of M D6'-branes, 
we would have a U (N) x U (M) gauge theory with the fermions supported at the inter- 
section points transforming as bifundamentals, and the above four-fermi interactions. 
To get an 577(5) theory with fields transforming in the 10, we should consider the case 
that N = M = 5 and the two stacks are exchanged by an orientifolding operation QR 
and intersect on the orientifold plane (f2 is the worldsheet parity; R acts by complex 
conjugating all three complex coordinates: Zi — > Zi, i — 1,2,3). Moreover, we pick 
an orientifolding operation that projects onto fermions whose wavefunctions A, and 
Xj transform as antisymmetric, rather than symmetric, tensors of SU(5). Tree level 
amplitudes in the orientifold theory are obtained by computing tree level amplitudes 
on the covering space for states that are invariant under the orientifolding projection, 
and then dividing by 2. The 10 2 10 2 interaction in an orientifold theory is thus derived 
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from 

4 



A A (k u k 2 , k 3 , k A ) =i(irg s a')I(6 l , 6 2 , 6» 3 )(27r) VQ^ ki)ui^u 2 u^^u A 



i=l 



,~ ~ ~ ~ x (33) 
Tr ( A1A2A3A4 + A1A4A3A2J • 



Evaluation Of I 

Now let us discuss the numerical evaluation of the integral I. After writing 

1 =i + A. (34) 



x(l — x) X 1 — X ' 

noting that the two terms contribute equally, and setting x = e~* to evaluate the 
contribution of the first term, we can write I as 

I(0i, 02, 9 3 ) = 2 r dt n JM^0ij[F(0 h l-0f,l; e-*)F(0 i , 1-^; 1; 1-e^)]- 1 / 2 . (35) 

Let us try to evaluate the integral on the RHS numerically for the most symmetric 
choice of rotation angles: 0i = 2 = #3 = 2/3. It turns out that for sufficiently large t, 
it is hard to maintain numerical precision in evaluating the hypergeometric functions. 
To deal with this problem, we will break up the integral into the range from to T 
and from T to 00. In the first region we use Mathematica to evaluate it numerically, 
while in the second we may use the asymptotics (20) to replace the integral by 

2tt 3 / 2 / dt(t + 3 In 3) " 3 / 2 = 4tt 3 / 2 (T + 3 In 3) ~ x l 2 , (36) 

JT 

where we used In 5(2/3) = 3 In 3. In practice, the sum of the two is insensitive to T in 
a certain range; we have checked this for T between 15 and 25. We find 

7(2/3, 2/3, 2/3) « 11.52 . (37) 

It is interesting that a large number 11.52 emerges from an explicit string calculation. 

Now cosider a generalization to 0\ = 2 = and 63 = 2 — 26. For the purpose of 
numerical evaluation, we approximate 1(6, 0,2 — 26) by 

T [F(2 - 26,26 - 1; 1; e' t )F(2 - 26, 26 - 1; 1; 1 - e^)]" 1 / 2 



2sm(7v6)^J-sm(27r6) dt 



F(6, l-6;l; e~ t )F(6, 1 - 6; 1; 1 - e"*) 

(38) 

+2- 3/2 * 



> T (t + \n5(6))Jt + ln5(2-26) 
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It is interesting to ask how 7 behaves when one of the rotation angles becomes small. 
This can be determined from (38) by letting 6 approach 1/2 from above, so that 
#3 = 2 — 29 corresponds to rotation by close to tt. It turns out that 7 decreases rather 
slowly from its maximum value at — 2/3. For example, 7(0.55,0.55,0.9) ~ 9.89 
and 7(0.51,0.51,0.98) ~ 6.69. Even when all three rotation angles become effectively 
small, 7 does not fall off very rapidly. Using (38), we find 7(0.9,0.9,0.2) 7.675; 
7(0.95,0.95,0.1) « 5.505; 7(0.99,0.99,0.02) « 2.47. We conclude that for a broad 
range of angles the quantity 7 lies in the range 7 — 11.5. 



4 Comparison To Four- Dimensional GUT's 

According to [35], eqns. (12.1.10b) and (13.3.24), the gravitational action for Type IIA 
superstrings is 

^j^^Ge^R. (39) 

The string coupling constant is g s = e*. After reduction to four dimensions on a 
six-manifold X of volume V 6 , the gravitational action in four dimensions is 

As the coefficient of 7? in four-dimensional General Relativity is conventionally (167rG n)^ 1 , 
with Gn Newton's constant, we have 

2 4 8V 6 G N 
g ° a (2tt) 6 1 ' 

The gauge coupling go6 of gauge fields on a 7J>6-brane is defined by saying that the 
effective action for the gauge fields is 



-V / <fx^Tr FijF ij , (42) 

where F^- is the Yang-Mills field strength and Tr is the trace in the fundamental 
representation of U(N). If we take the 7J>6-brane worldvolume to be R 4 x Q, where Q 
is a compact three-manifold of volume Vq, then the action in four dimensions becomes 



V, 



Q 



^9 



, J d 4 xTr FijF ij . (43) 



DG 



However, before comparing to four-dimensional GUT's, we must recall that it is conven- 
tional to expand the gauge fields as Ai = J2 a QaA°;, where Tr Q a Qb = \$ah- Similarly, 
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one conventionally expands F^ = J2 a FijQa- If this is done, the action becomes 

Vn 



, , d*xJ2F«F^. (44) 

The GUT action is conventionally written 

[d 4 xJ2F^ a , (45) 

where Qgut is the GUT coupling. Hence, we have g GUT = 2g 2 m /VQ. Since g 2 m = 
(2n) 4: g s a' 3 ^ 2 according to eqn. (13.3.25) of [35], we have 

2(2ir)%a' 3 / 2 

9 GUT = 77 • (46) 

V Q 

The volume Vq enters in the running of SU (3) x SU (2) x U (1) standard model gauge 
couplings from very high energies down to the energies of accelerators. Roughly speak- 
ing, Vq 1 ^ 3 plays the role of M GUT , the mass scale of unification, in a four- dimensional 
GUT theory. To find the precise relation between them, one must compute the one- 
loop threshold corrections to gauge couplings. This was done in [12] for M-theory on 
a manifold of G2 holonomy. The one-loop threshold corrections are the same in Type 
IIA as in M-theory, since they come from Kaluza-Klein harmonics on Q that are the 
same in the two theories. So we can borrow the result of [12]. According to that result 
(see eqn. (3.30) of [12]), the precise relation between M GUT understood as the mass at 
which the low energy coupling constants appear to unify and Vq is 

ivi GUT 

where L(Q) is a certain topological invariant of Q (together with Wilson lines on Q that 
break the GUT symmetry to the standard model) that is known as the Reidemeister 
or Ray-Singer torsion. L(Q) depends on a model, but is readily computable in a given 
model. For example, in a model considered in [12] in which Q is a lens space S 3 /Z (? for 
some positive integer q (which must be prime to 5 if we want an SU(5) model), and 
the Wilson line on Q has eigenvalues exp(2iri5i/q) with Si = (2w, 2w, 2w, — 3w, — 3w) 
for some integer w prime to q, one has 

L(Q) = 4gsin 2 (5vrw/g). (48) 

For a slightly more general lens space, also described in [12], whose definition depends 
on another integer j, this is replaced by L(Q) = Aq\sin(5nw/q)sm(5njw/q)\. The 
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factor L(Q) can be important in determining the proton lifetime, as will become clear. 
For example, for the minimal choice that leads to the standard model gauge symmetry, 
which is q — 2, w — 1, we have L(Q) = 8, and the threshold correction will prove 
to enhance the proton decay rate significantly. It is possible for lens spaces to make 
L(Q) < 1, but only if q is extremely large. 
Using (48), we can express (46) in the form 

a GUT = (2irfL(Q)- l g s a ,V2 Ml UTl (49) 

where as usual a G ur = QgutI^ 11 - Alternatively, 

2 /3 _ <*GUT L (Q) 2 (t -n\ 

Ideally, we would like to compute all of the parameters in the string compactification 
in terms of the quantities G^ = M GUT , and a G uT, about which we have at 

least some experimental knowledge. The Planck mass is well-known {M Pi = 1.2 x 10 19 
GeV), but M GUT and oi G ut are somewhat model-dependent. The most commonly 
cited values based on extrapolation from low energy data are M G ut = 2 x 10 16 GeV, 
a G uT — -04. Unfortunately, the string theory really depends on four unknowns V 6 , Vq, 
g s , and a'. 

To parametrize our ignorance, we may introduce the dimensionless quantity 



A = ||, (51) 



which is of order one if X is fairly isotropic, and solve for the other stringy parameters 
in terms of Gn, M gut , a G uT, and A. For Vq, this has already been done in (47). 
Dividing (41) by (50), we get 

, 8\G N 

a = — 2 . (52) 

a GUT 

We can similarly solve for g s , 

9s = a ou T L{Q) 

2 9 / 2 (2vr) 3 M^ T < /2 A 3 / 2 

The factor about which we have the least intuition is A, which comes from a ratio of 
volumes. To try to get some intuition about the possible values of A, and also about 
whether the model makes sense, let us examine quantitatively the formula (53) for g s : 

a n i f aGUT Y f 2xl0 16 GeV \ 3 L(Q) 
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Our calculation does not make much sense if g s ^> 1, since then we should do the 
computation in M-theory (leading back to the discussion in section 5 of [12]) rather 
than in Type IIA superstring theory. If g s <C 1, our computation makes sense, but it 
is implausible for the vacuum to be stabilized after supersymmetry breaking in a way 
that would enable the world as we see it to exist. So the discussion makes most sense 
if g s is relatively close to 1. We note that happily (54) is compatible with having g s 
near 1 while the GUT parameters have their usual values, L(Q) = 8, and A is not too 
far from 1. 

In view of the preceding discussion, it is perhaps more useful to use (54) to solve 
for A, parametrizing our ignorance via the unknown value of g s , which we expect to be 
near 1: 

x _ 8/3 Ml, L(Q) 2/3 
~ aGUT M* UT s(2nygr (M) 

Now (52) becomes 

(27r) 2 ^ 2/3 M ( 

Finally, the amplitude for proton decay is, from (33), 



a' = "out^j (56) 



A st — Tra g s l{Ui, V 2 , V 3 ) — 2 , 

AttM gut 

where / is the integral discussed in the last section. It is interesting to compare this 
result with the M-theory estimate of [12]: 

af UT L{Q)^ 

A M 772 • ( 58 ) 

ivi GUT 

We find the same scaling with a GUT) L(Q) and M GUT as in M-theory. If we keep these 
parameters and the rotation angles #j fixed, then we may study the amplitude as a 
function of g s . String theory indicates that this function behaves as g l J z for small g s ; 
M-theory tells us that it approaches a constant for large g s . Strictly speaking, (57) is 
reliable for small g s , but we will take it as a rough approximation for g s of order 1. 

Analogous Field Theory Amplitude 

Let us compare this to the analogous four-fermion amplitude in four-dimensional 
GUT's. For simplicity, we will assume that (as in SU (5) models) all superheavy gauge 
bosons have the same mass M X - M x is comparable to the unification scale M GUT 
inferred from the running of the low energy gauge couplings, but differs from it, in 
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general, by model-dependent factors. Exchange of superheavy gauge bosons of mass 
X gives an amplitude 

^E(^I^|A 2 )(A 3 |^|A 4 ), (59) 

where we have labeled initial and final fermion states as A 1 , A 2 , A 3 , and A±. The sum 
runs over the superheavy gauge bosons, but since the standard model gauge bosons 
contribute baryon-number conserving interactions anyway and taking the sum over all 
generators of the gauge group will lead to a simpler formula, we will do this. (Standard 
model gauge boson masses are near zero, not near Mx, so (59) is not a good approx- 
imation to their contribution.) Note that, ignoring worldsheet instantons, the string 
theory four-fermion amplitude is S'?7(5)-invariant, since it comes from a local contribu- 
tion that does not see the symmetry breaking by Wilson lines, while the four-fermion 
operator in GUT's is of course not S*?7(5)-invariant. 

We first consider the case that the fermions transform as 5's of SU(5), though 
the resulting 5 2 5 2 amplitude is actually baryon number conserving. To get baryon 
nonconservation, we need to incorporate 10's as well, as we will do presently. Also, to 
get a better match to the local string theory construction with Chan-Paton factors, we 
take the gauge group to be U(5) instead of SU(5); the extra U(l) gauge field does not 
violate baryon number anyway. (In string theory, the local construction at a particular 
brane intersection point has U(5) gauge symmetry, but globally in realistic models the 
extra U(l) is Higgsed by absorbing an RR mode.) 

Let us work out the group theory part of the matrix element in eqn. (59). For this, 
we introduce column and row vectors and ctj for states transforming in the 5 or 5; 
the group theory factor of the current-current matrix element becomes 

J2aiT a a 2 a 3 T a a 4 , (60) 

a 

where the T a are traceless 5x5 matrices that generate U(5). As we noted earlier, the 
T a are conventionally normalized so that Tr T a Tb = \5 a b- It is straightforward to show 
that 

E(aiT a a 2 ) (a 3 T a a 4 ) = -(a 1 a 4 )(a 3 a 2 ). (61) 

a * 

The 5 2 5 2 interaction in GUT's is hence 

9g 

T^pP [(aia A )(a 3 a 2 ) + (aia 2 )(a 3 a 4 )] , (62) 

where we sum over the two channels in which the gauge boson can be exchanged. 

Now we move on to the 10 2 10 2 interaction. In order to make clear the origin of a 
certain factor of 2, we will introduce the 10 by a field theory version of the orientifold 
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that we used in string theory. The gauge group is £7(5) x £7(5), with gauge fields A, 
A' and the fermions consist of fields ip, ip' transforming as (5, 1) + (1, 5) plus a field S 
transforming as (5, 5). The Lagrangian is invariant under an "orientifolding" symmetry 
that exchanges A with (A') T , where (A') T is the transpose (or equivalently, as A' is 
hermitian, the complex conjugate) of A', breaking £7(5) x £7(5) to a diagonal £7(5). © 
exchanges ip and ip', leaving a 5 of the unbroken £7(5), and we take it to act on S with 
a suitable sign such that the invariant modes in S transform as the 10 of the diagonal 
£7(5), corresponding to an antisymmetric 5x5 matrix S %J . 

The kinetic energy is taken to be a general 9-invariant expression: 

/ = —j— Tr(F(A) 2 + F(A'f ) + ^7 • Dip + ^7 • Dip' + Si') ■ DS. (63) 



4<?, 



hi 



where sums over all indices of ip, ip', and S are understood. Let us work out the S 2 S 2 
interaction prior to orientifolding. It comes from exchanges of gauge bosons in the two 
£7 (5) groups. In working out the contribution from exchange of a gauge boson in either 
£7 (5) factor, the indices of S that transform under the other £7 (5) are spectators. The 
amplitude due to exchange of either £7(5) is thus just like (62), except that when we 
include the spectator index, the wavefunctions all are 5x5 matrices \ and Xj rather 
than row and column vectors, and the pairing of row and column vectors is replaced 
by a trace. The amplitude due to exchange of either £7 (5) is hence 

9GU 2 T ^17^2^37/^4 Tr(AiA 2 A 3 A 4 + AiA 4 A 3 A 2 ). (64) 



2M GUT 



The total amplitude due to exchange of a gauge boson in one or the other £7 (5) is hence 

g 2 .... 

• G 2 T ^i7^^2%7m m 4 Tr(AiA 2 A 3 A 4 + AiA 4 A 3 A 2 ). (65) 



M GUT 



Orientifolding is carried out by restricting to G-invariant states and dividing the 
action by 2. As a result, we can drop A' and ip', but we must divide the S kinetic 
energy by 2. The resulting kinetic energy is then 

V = AgJ^ F(A)2 + ^ 7 ' D ^ + l^iW ■ DS ij . (66) 

The S kinetic energy has been divided by 2, but is still canonically normalized in the 
following sense. S is an antisymmetric matrix, so for example has a 2 x 2 block looking 
like 

(67) 



e 
-e+ 
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The kinetic energy is canonically normalized for fields like e + appearing above the 
diagonal in S. 

As for the 10 2 10 2 interaction in the orientifold theory, we can borrow it from (65). 
We merely have to take the wavefunctions Aj and Xj to be invariant under O, which 
means that they are 5x5 antisymmetric matrices. Also, we have to divide by 2, just as 
we divided the O-invariant part of the classical action by 2 to get (66). So the 10 2 10 2 
interaction is finally 

r wi7^2«37^4 Tr(AiA 2 A 3 A 4 + AiA 4 A 3 A 2 ). (68) 



2Mq UT 



Of course, we could have derived this directly starting with the kinetic energy (66). 
Comparison 

To compare string theory to four- dimensional gauge theory, we therefore need merely 
compare the string theory amplitude ng s a'I from (33) to the field theory amplitude 
9gut/2M x = 2nacuT / M x from (68). The comparison is thus 



7cg s al ^2ttj^. (69) 



We call the left hand side the string amplitude A st , and the right hand side the corre- 
sponding GUT amplitude Agut- Using (57), we find the ratio to be 



A st L{QflH g]/z M 



~ X (70) 



Agut " 8vr 2 af UT M GUT 

1 /3 

For other parameters fixed, this ratio scales as l/o; Gl/T , which is the same as in the M- 
theory calculation of [12]. The negative power of cxgut reflects the fact that in this type 
of model, proton decay is a stringy effect, rather than being dominated by exchange 
of a single Kaluza-Klein mode. Exploring a possible enhancement due to this factor 
was the motivation for the present paper, but the factor of 87r 2 in the denominator is 
clearly unfavorable. We can rewrite (70) as 

A« „,„_ / .04 X 1 / 3 Ml 



0m7L(Q)*"lgW (-^-) (71) 



Agut s \&gut/ M gut 

Happily, g s , of which we can only say that we expect it to be more or less close to 1, 
is here raised to a relatively small power, helping to reduce the uncertainty. 

Using the most commonly cited values Mqut = 2 x 10 16 GeV, acur — -04, a recent 
evaluation of the proton lifetime in four-dimensional SUSY SU (5) due to gauge boson 
exchange gave the value 1.6 x 10 36 years if Mx = Mqut [5] (we have taken eqn. (17) 
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of this paper with M x = 2 x 10 16 GeV). If M x or a GUT is changed, the proton lifetime 
scales as: 



Using modified values of (Xgut and Mqut is natural in the present discussion if doublet- 
triplet splitting is accomplished with discrete symmetries as in [8], since this necessi- 
tates extra matter fields (in complete SU(5) multiplets) surviving to energies below 
the GUT scale. For example, in [36], in models with such additional multiplets, it was 
found plausible to have values such as a GUT = -2, M GUT = 8 x 10 16 GeV. The net 
effect of these changes is to roughly double the proton lifetime. 

To obtain the proton lifetime r p si in this class of string theories from (72), we multiply 
by (Agut / A st ) 2 and replace the prefactor 1.6 by 2 for a reason explained below. Thus: 



r p „ t = 2 X ^ (£-f (^^vf ■ (73) 



A proton lifetime of order 10 36 years due to gauge boson exchange is considered 
unobservably small for the foreseeable future; the present experimental bound on p — > 
7r°e + is about 4.4 x 10 33 years (for example, see [37]), while a next generation detector 
(described for example in [38]) may reach a limit close to 10 35 years. However, it is 
clear that a relatively modest enhancement in the proton decay amplitude might save 
the day. Looking at (71), we see a few large and small factors that tend to cancel. This 
formula contains the small overall factor 0.037, along with the factor /, which we have 
found to be roughly 7 — 11 for plausible values of the angles, and the threshold factor 
L(Q) 2 / 3 , where L(Q) need not be large but is in fact equal to 8 for the minimal lens 
space with fundamental group Z2, and a little larger for many of the lens spaces with 
other small fundamental groups. Combining these factors, the best we can say is that 
in a model based on intersecting D6-branes, rather as in four-dimensional GUT's, the 
proton lifetime due to dimension six operators is likely to be close to 10 36 years. 

Note that one major source of uncertainty in GUT's is absent here: the proton life- 
time directly involves the scale M GUT that can be probed using low energy data, rather 
than the heavy gauge boson mass M x whose relation to M GUT is model-dependent. 

What we have evaluated in this paper is a local, stringy contribution to proton 
decay in a certain class of models based on intersecting Type II D6-branes. If the 
compactification volume is even slightly larger than the string scale, so that worldsheet 
instanton effects are at least slightly suppressed, what we have computed is likely to 
be the dominant contribution (from dimension six operators) in this class of models. 
The result we have obtained has an anomalous power of oi G ut because it is a short 
distance stringy effect. As in [12], this contribution yields only a 10 2 10 2 interaction, 




(72) 
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which contributes to p — > 7r°e^. In four-dimensional GUT's, there is also a 1010 5 5 
contribution, leading to p — > 7r°e^. With the assumptions made in [5], the ratio of 
p -> 7r°e^ to the total is y = 1/(1 + (1 + |Kd| 2 ) 2 ), where V ud = 1 is a CKM matrix 
element. Hence in comparing the stringy proton decay rate to that in field theory, we 
should include a factor of 1/(1 — y) in the proton lifetime in the string model to account 
for the missing p — > vr°ej. We included this factor in taking the prefactor in (73) to 
be 2. 

The last paragraph has been formulated loosely; with different assumptions about 
the flavor structure at the GUT scale, the tc° could be a K° and the e + could be a /i + 
(in which case the lepton polarization in the final state would be measurable). In either 
the GUT theory or the string model, the proton lifetime could be increased further by 
unfavorable assumptions about flavor structure (mixing with the third generation, for 
instance, or mixing with new GUT-scale fermions). At any rate, the assertion of this 
class of models that proton decay is caused mainly by 10 2 10 interactions is testable 
in principle. 
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